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We study the curvaton scenario using gauge-invariant second order perturbation theory and solv- 
ing the governing equations numerically. Focusing on large scales we calculate the non-linearity 
parameter /nl in the two-fluid curvaton model and compare our results with previous analytical 
studies employing the sudden decay approximation. We find good agreement of the two approaches 
for large curvaton energy densities at curvaton decay, fio-dcc, but significant differences of up to 10% 
for small f^crdcc. 
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I. INTRODUCTION 

^^ ' The third year WMAP data release [l| has confirmed beautifully the cosmological standard model of structure 

ZTj formation: during inflation fluctuations in the scalar fields are stretched to super-horizon scales and later on source 

the Cosmic Microwave Background (CMB) anisotropics and the large scale structure. In the standard inflation 

j^ models the scalar field responsible for the accelerated expansion of the universe, the inflaton, also provides these 
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fluctuations 13- Recently a related but different scenario has become become popular: the curvaton paradigm 



j-^ , nuctuations la- tiecentiy a reiatea C 

0, 0, IS la H la S Ha HH uM- Here the fluctuations are not generated by the inflaton, but by a different scalar 
Cn ' field, the curvaton. 



A powerful tool to differentiate between different models of the early universe is second order perturbation theory 



A poweriul tool to dmerentiate between ainerent models oi tne early universe ] 
f^'. ^M,M,M,M,M,M,^M,^^^M,M,^M,^M,M,^^ it allows the accurate calculation 
^ . of higher order statistics such as the primordial bispectrum and the non-linearity parameter /nl- Instead of using 
[^ ' cosmological perturbation theory at second order, the AA^-formalism _34, _35^, _36j _37^, _38j _39J has proved useful to 
OO . study non-gaussianity and calculate /nl- 

C^ ■ Bartolo et al. ^(f\ studied the curvaton model at second order using cosmological perturbation theory, and Lyth 

"nI^ ] and Rodriguez used the AiV-formalism to calculate the non-linearity parameter /nl in this scenario. However, both 
studies used the sudden decay approximation. Here we go beyond sudden decay, using second order gauge-invariant 
perturbation theory and solve the ensuing equations numerically. 

We consider scalar perturbations up to and including second order and assume a flat Friedmann- Robertson- Walker 
i-Q ■ (FRW) background spacetime. We work on large scales (compared to the horizon size), which allows us to neglect 
f^. gradient terms. 

Q ■ The outline of the paper is as follows. In the next section we give the governing equations up to second order 

H [ and deflne the gauge-invariant variables we are using. We specify in Section IIIII the two-fluid curvaton model we 
c/3 ■ are studying and apply the equations of Section ^] In Section IIVI we take a small detour from perturbation theory 
and compare our perturbative approach to the AN formalism. After deflning the non-linearity parameter /nl we 
present numerical solutions in Section and compare our numerical results to the sudden decay approximation. The 
governing equations without any gauge restrictions are given in the appendix. 



II. GOVERNING EQUATIONS 

In this section we give the governing equations for a system of multiple interacting fluids on large scales, allowing 
for scalar perturbations up to second order, following closely the treatment of Refs. |23 and |4(]|. 
The covariant Einstein equations are given by^ 

G^, = 8^G T^, , (2.1) 



^ Notation; Greek indices, fj,, u, A, run from 0, ... 3, while lower case Latin indices, i,j, k, run from 1, ... 3. Greek indices from the beginning 
of the alphabet, a, fi, 7 will be used to denote different fluids. 



where G^^ is the Einstein tensor, T^^ is the total energy-momentum tensor, and G is Newton's constant. Through 
the Bianchi identities, the field equations H2.1|l imply the local conservation of the total energy and momentum. 



^^j^f^i^ = 0. 



(2.2) 



In the multiple fluid case the total energy-momentum tensor is the sum of the energy-momentum tensors of the 
individual fluids 
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(2.3) 



For each fluid we deflne the local energy-momentum transfer 4-vector Q'f^) through the relation 



(2.4) 



where energy-momentum is locally conserved for Q\^) = 0, i.e. only for non-interacting fluids. Equations l|2.2|) 
and H2.4|l imply the constraint 
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We split scalar perturbations into background, first, and second order quantities according to, 

using here the total energy density as an example. 
The line element on large scales is given by 



(2.6) 



ds^ 



I + 2 [ 4)1 + -(f>2 



dr + a^ 



1 - 2 ( ^1 - -V2 



Sijdx^dx-' 



(2.7) 



where a = a(t) is the scale factor, 4>i and 02 are the lapse functions at first and second order, respectively, and ipi 
and ip2 the curvature perturbations. 

Following Refs. p(ll4ll| we split the energy-momentum transfer 4-vector using the total fluid velocity u^ as 
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where Qia) is the energy transfer rate and fi^ s the momentum transfer rate, subject to the condition u^f^ ^ = 0. 



On large scales the only non-zero component of the 4- velocity is 

1 , 1 



uo = 



l+(Pl^ -(P2 ~ 2' 



(2.9) 
(2.10) 



subject to the constraint m^u^' = — 1. We then find the only non-zero component of the energy transfer 4-vector on 
large scales to be 



Qo{a) = -Qoa ( 1 + 01 + -02 " -^(pj ] - SQla (1 + 0l) - 
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(2.11) 



where Qoaj ^Qia, and SQ2a are the energy transfer to the a-fluid in the background, at first and at second order, 
respectively. 



A. Background 

Energy conservation for the a-fluid in the background is given from Eq. 12.41) as 

POa ~ — 3iJ {poa + Poa) + Qoa , 



(2.12) 



where H = a/ a is the Hubble parameter, and poa a-nd Pqq are the energy density and the pressure, respectively, of 
the a-fluid. Total energy conservation is then given by summing over the individual fluids and using Eq. (|2.5|l as 

po = -3H (po + Po) , (2.13) 

where Po = LaPoa- 

The Friedmann constraint is given from the — component Eq. H2.1|l as 

i^^ = ^Po. (2.14) 

B. First order perturbations 

We now give the governing equations on large scales at first order in the perturbations in the flat gauge, denoting 
quantities evaluated in this gauge by a "tilde" . The governing equations in an arbitrary gauge are given in appendix 

Em 

The energy conservation equation for the a-fluid at first order is given from Eq. (|2.4|) on large scales as 

(57]^ + 3i/ (^7^ + <5F|^) - Qoa^T - (5g;^ - . (2.15) 

The total energy density perturbation is related to the individual fluid densities, and similarly for the pressure 
perturbations, by 

j;^ = ^jTTo, SK^Y.^-" (2.16) 

a a 

and we get, using the constraint Eq. I|2.5|l . from Eq. (|2.15() the evolution equation for the total energy density 
perturbation 

Jpi+iH ijf^i + SK) = . (2.17) 

The — Einstein equation on flat slices is, using Eq. H2.1|l and the background Friedmann constraint H2.14|l . given 
by 



2 Po 
The curvature perturbation on uniform a fluid energy density hyper-surfaces at first order is given by |42 



(2.18) 



Ci„ = -ff^. (2.19) 

The curvature perturbation on uniform total energy density hyper-surfaces at first order is given by 

Ci=-?-^, (2.20) 

Po 

and related to the curvature perturbation on uniform a fiuid slices by 

Ci = E— Ci- (2-21) 

^ PO 

As in the background we introduce new variables, the normalised energy densities at first order, 

m^^ = ^ , (2.22) 

Po 

which allow us in combination with choosing the number of e-foldings as a time variable to write the governing 
equations in the following sections in a particularly compact form. 

In terms of the new variables the curvature perturbation on uniform a fluid energy density hyper-surfaces, given 
in Eq. H2.19(l . is simply 

Cia = -H-^5nia, ■ (2.23) 

POa 



C. Second order 



We now give the governing equations on large scales at second order in the perturbations in the flat gauge, denoting 
quantities evaluated in this gauge by a "tilde" . The governing equations in an arbitrary gauge are given in appendix 

EH 

The energy conservation equation for the a-fluid at second order is given from Eq. H2.4|l on large scales by 

-VTa + iH (Jf^„ + sKl) - Qoa (fi ~ Jx^ - SChc - 20i5'Q^ - . (2.24) 

Using Eq. 1)2. 5|l the evolution equation for the total energy density is 

'5p2 + iH{ip2 + 5P2)^Q. (2.25) 

where the total density and pressure perturbations are given in terms of the individual fluid ones by 

^2 = Y,^'^' S^ = ^^. (2.26) 



by 



The — Einstein equation on flat slices is, using Eq. H2.1|l and the background Friedmann constraint H2.14(l . given 



\5p2_ 
'2 po 



(2.27) 



The curvature perturbation at second order in terms of uniform a-density perturbations on flat slices is given by 
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(2.28) 



where c^ = Poa/poa is the adiabatic sound speed of the a-fluid. 

The curvature perturbation at second order in terms of the total density perturbations on flat slices is given by |2E 



H 



H 



H' 



C2 = — —Sp2 + 2^6pi6pi + ^2" (5 + ^cj 5pi , 
Po Pt; Pa 

where Cg = Pq/po is the total adiabatic sound speed related to the individual speeds c^ by 

2 Y^ POa 2 
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(2.29) 



(2.30) 



As at first order we introduce new variables, the normalised energy densities at second order, allowing us to rewrite 
the governing equations in the following sections in a particularly compact form. 
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(2.31) 



The curvature perturbation (^2, defined above in Eq. H2.29|l . is related to the curvature perturbation employed in 
the AA^ formahsm (see also |43|), which we denote by C2SB, by |2Ml37j 



C2SB — C2 — 2^1 



(2.32) 



It was originally introduced by Salopek and Bond |3J| and employed by Maldacena in studies of non-gaussianity in 
Ref. El. 




FIG. 1: Evolution of the normalised background curvaton density, Qotr, and the normalised decay rate g, as a function of the 
number of e-foldings, starting with initial density and decay rate fiotr ~ 10^^ and F/H — 10"'^, corresponding to pin = 0.32 on 
the left panel, and on the right panel with F/H = 10~^, corresponding to pin = 10. 



III. THE MODEL 

In this section we specify the curvaton model and applythe governing equations given in the previous section order 
by order. We model the curvaton as a pressureless fluid [3j, and hence our system will be governed by the equations 
of state 



J'^a , -Tj — q ^7 ' 



(3.1) 



where the subscripts "cr" and "7" denote the curvaton and the radiation fluid, respectively. The decay of the curvaton 
is described by a fixed decay rate, F — const, 



Qa = -^Pa , Q-1 = Tp, 



7 -^ Pct j 



where we used Eq. H2.5|l . 



(3.2) 



A. Background 

The background evolution equations arc from Eq. H2.12|l and using Eqs. H3.1|l and H3.2|l given by 

Pa = -'iHpQa - EpoCT , (3.3) 

P07 = -^Hpo^ + Tpoa ■ (3.4) 

We now change to a new set of variables. First we introduce normalised energy densities in the background. 



O - ^°" 

"Oq = 



and define the reduced decay rate as 



V + H 



(3.5) 



(3.6) 



We change the time coordinate from coordinate time t to the number of e-foldings iV = Ina, that is ^ = ^'IN- "^^^ 
normalised radiation energy density is then simply given from the Friedmann equation, l|2.14(l , as VIq^ = 1 — IIqit and 
we get the system of background evolution equations in terms of these new variables 



i^n^ — f^nrr 1 — ^nn — 



'Oct 



1-5 



g' = -(4-r!o.)(l-5)5. 



(3.7) 
(3.8) 



Solutions for the system H3.7|l and (|3.8() are given in Fig. ^ for two different initial conditions, fJocr = 10 ■^ and 
T/H = IQ-^ and n^a = lO^^ and T/H = lO^^. 
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FIG. 2: Evolution of the total curvature perturbation, C,\, and the normalised density perturbations at first order as a function 
of the number of e-foldings, starting with i^io- = 1 and initial density and decay rate noo- = lO"'^ and V jH = 10~^, corresponding 
to Pin = 0.32. 
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FIG. 3: Same as Fig. |21but with V/H = 10 ® initially, corresponding to pin — 10. 



It was shovifn in Ref. j7| that the solutions of the system (|3.7() and (|3.8() depend only on a single parameter since 
we can write ^oa = ^^i3)^ and for t; ^ 1 we can solve the system explicitly, which gives J7o<j oc y^. We therefore 
define the parameter 



Pm 






the subscript "in" denoting the initial conditions. For the initial conditions r^oo- = 10 ^ and T/H 
r2o(T = 10~^ and T/H = 10"^ the parameter pin takes the values 0.32 and 10, respectively. 



(3.9) 



10^3 and 



B. First order 



The perturbed energy transfer rates are given from Eq. 1)3. 2|l at first order as 

SQla = -T6pia , 5Qi~f = TSpia ■ 



(3.10) 



The evolution equations at first order are from Eq. (|2.15() and Eqs. I|3.1|l . and using l|3.10(l in terms of the normalised 
energy densities defined in Eq. H2.22|l given by 
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(3.11) 
(3.12) 
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FIG. 4: Evolution of the total curvature perturbation, C,i, and the normalised density perturbations at second order as a 
function of the number of e-foldings, with initially ^lo-.in = 1 and density and decay rate ^0^ = 10~^ and V jH — 10"'^, 
corresponding to pin = 0.32. 



The curvature perturbations at first order on uniform curvaton and radiation density hypersurfaces are from 
Eq. H2.19|l given by 

Clcr = —, -^TTT^ — ^^\c 
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l7 
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4(1- g) Qqj - gQo„ 
The curvature perturbation on uniform total density slices in terms of the new variables is given by 
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Afl, 



O7 



(3.13) 
(3.14) 

(3.15) 



Solutions for the equation system H3.11|l and (|3.12|l are given in Figs. [3 and |21 for two different sets of initial 
conditions, floa = 10~^ and T/H — 10~^ and Uq^ — 10^^ and T/H = 10~^, corresponding to pin = 0.32 and pin — 10, 
respectively. Using Eq. I|3.15|) we also plot the evolution of Ci ■ For the perturbations we use the initial conditions 

ClCTjn = 1 , Cl7,in=0, (3.16) 

which can be easily translated in initial conditions for Silia using Eq. H3.13|l . and facilitates comparison with Ref. J2f . 
Note that the values for Sflia and 6fli^ can exceed 1, as can be seen in Fig.|31 This doesn't indicate the "breakdown 
of perturbation theory" or anything dramatic like it, but is merely an artifact of normalising the density perturbations 
by the total background density po, which can itself be small. As in the background, the normalised energy densities 
together with the choice of time coordinate give a particularly neat system of governing equations. 



Second order 



The perturbed energy transfer rates are given from Eq. H3.2|l at second order as 

SQ2a = -^Sp2a , SQ2'y = T6p2a ■ (3-17) 

We then find evolution equations at second order from Eq. H2.24|l and using Eqs. H3.1|) and (|3.17|l in terms of the 
normalised energy densities defined above in Eqs. H2.22II and H2.31|l to be 
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FIG. 5: Same as Fig. |l|but with F/H = 10 ^ initially, corresponding to pin ~ 10. 



The adiabatic sound speed in a multi-fluid system is given above in Eq. 1)2. [^0(1 and we find for the two-fluid curvaton 
model 



1 4(1-5) ^oy - 5^ 
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The curvature perturbation on uniform total density hypersurfaces at second order in terms of the normalised 
quantities is 
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(3.21) 



where Cg is given above in Eq. 1)3.20(1 . 

The system of equations ((3.1811 and ((3.19(1 is readily integrated using a standard fourth order Runge-Kutta solver 
p4l |. We give the solutions for this system of equations for the two different sets of initial conditions, floa — 10~^ and 
r/i/ = 10-3 and Qoa = 10"^ and T/H = 10-^ corresponding to pin = 0.32 and pi„ = 10, in Figs. HandEl Using 
Eq. 1(3.21(1 we also plot the evolution of ^2. The initial conditions for the second order perturbations are chosen as 
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Note that the values for Sil2a and Sfl2j can exceed 1, as can be seen in Fig. O As at first order, this is merely an 
artifact of using normalised energy densities. We do however see a new effect: at second order the energy densities 
Sp2a can and do become negative, as can be seen clearly in Figs.0]and[Hl This is a "real" effect and not a normalisation 
artifact since po ^ always. However, the total energy density p as given by summing over all the terms in the power 
series expansion Eq. 1(2.6(1 . again stays positive definite. 

IV. RELATING THE PERTURBATIVE TREATMENT TO THE AN FORMALISM 



The AA^ formalism |35| provides a simple tool to calculate the curvature perturbation on large scales at all orders 
in the perturbations on scales larger than the horizon |3a, |33, ISg ■ The main simplification compared to cosmological 
perturbation theory stems from the fact that we only need the background evolution equations, and not the full 
governing equations at all orders of interest. However, if there is no analytic solution the numerics necessary to get a 
result turn out to be quite involved as can be seen below. Nevertheless, we shall outline the calculation in the following. 

The AA^ formalism relates the curvature perturbation on uniform density hypersurfaces to the perturbation in the 
number of e-foldings from the uniform density to the flat slicing. 



To get the number of e-foldings A^ we use Eq. ((3.8() to get dN in terms of dg, and integrate, 
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The curvature perturbation in the AA'^ formahsm is then given from Eq. H4.1|) by expanding A'^ in a Taylor series, 
which leads in the curvaton case to \12L l38l | 
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where the partial differentials are 
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Note, that to make contact with first and second order perturbation theory, C and (Jrio-in have to be expanded up 
to second order, where in this case the second order curvature perturbation corresponds to C2SB7 defined above in 
Eq. H2.32|l . related to C2 as specified in Eq. (|2.29() . 

Although in principle we can evaluate the integrals in Eqs. (|4.4I) and (|4.5|) numerically and then differentiate them 
with respect to the initial conditions to get the value of C, this is (arguably) more difficult than solving a set of 
coupled differential equations. We therefore don't use the AA^ formalism in the following sections, and solve instead 
the system of differential equations presented in Section ITTll 

However, the AA^ formalism is used in Ref. |45| in another numerical study of the curvaton scenario. The results 
are similar to the ones presented in this paper, but the computing time required in the AA'^ case is increased by factor 
of roughly ^ 100 compared to solving the system of differential equations presented in Section ITTTI 



V. THE NON-LINEARITY PARAMETER /nl: RESULTS AND DISCUSSION 



In this section we give the non-linearity parameter /nl calculated numerically using the governing equations at 
second order of Sections |n] and IIIII and compare it to previous numerical first order results and analytical sudden 
decay estimates. 



The non-gaussianity parameter /nl is defined as 



C = Cg + g/NL 



^S Sg 



(5.1) 



where (^ is the curvature perturbation at all orders, Cg the gaussian part of C, and the "bar" denotes the spatial 
average. There has been some confusion in the literature as to the sign of /nl, which becomes relevant if the result 
is compared with observations. The sign convention chosen here coincides with the one used originally by Komatsu 
and Spergel 46] and adopted by most observational studies, and corrects the sign error introduced in Ref. 6] and 
carried through in much subsequent work ^0, 24] ^. 



We now relate the curvaton field fluctuations to the curvaton fluid energy density. The energy density in the 
curvaton field can be approximated by 
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(5.2) 



where m is the curvaton mass and a is the amplitude of the curvaton field. Expanding the curvaton amplitude to 
first order, a = ctq + <5cri, we get from Eq. (j5.2|l . 
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^ Note that Eq. (36) of Ref. I3I, corresponding to Eqs. 15.81 and i5.12i here, has the correct sign for /n 
the derivation in Ref. la. 



however there is a sign error in 
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FIG. 6: The final values of the normalised curvature perturbations at first and second order, (^i,in and C2,in, for an initial value 
CiCT.in = 1, versus the background curvaton energy density at decay, Qcrdcc, in the left panel and versus pin in the right panel. 



Note, that including the quadratic term ((5cri/cro)^ hi the first order energy density and setting the second order energy 
density perturbation to zero is just convention, following Ref. |g]. In Ref. [23| this term is included in the second order 
energy density perturbation of the curvaton fluid. This choice doesn't effect the final results. 

A. Sudden decay 

In the sudden decay one assumes that the curvaton doesn't decay into radiation until a time idee, when all of the 
curvaton energy density decays suddenly; the normalised energy density of the curvaton at decay is denoted Oo-dee- 
The sudden decay approximation has been widely used in the literature to study the curvaton scenario without 
having to resort to numerical calculations, see e.g. Refs.0, 0, '20, 38]. In order to be able to compare the sudden 
decay approximation to the numerical calculation we now give a prescription to calculate rjg-dec- 



The normalised background energy density can be approximated by [21, llll Il2| 
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6M^J^ 



(5.6) 



where Mpj — (SttG) ^ . Assuming that the evolution of the background curvaton amplitude from the initial time up 
to curvaton decay is negligible and using that initially _ffi„ = m, and at curvaton decay i?doe — T we can use Eq. H5.6|l 
to relate the parameter pin, defined in Eq. (|3.9|l . to the background energy density of the curvaton at decay. We 
therefore define 
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as energy density of the curvaton in the sudden decay approximation. The agreement fla-dcc defined in Eq. H5.7|l with 
^crdee uscd in Rcf. 7] is quite good, and we use the definition H5.7|l in the following to compare our numerical results 
with the sudden decay approximation. 



We now briefiy review the results of previous analytical treatments using the sudden decay approximation to 
calculate the non-linearity parameter /nl in the curvaton scenario. 

The non-linearity parameter in the sudden decay approximation using first order perturbation theory, however 
using the definition of the first order energy density perturbation quadratic in the curvaton fluctuations Eq. H5.4|l . is 
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Usiiig second order perturbation theory the non-linearity parameter in the sudden decay approximation was found 
to be m,E^ 
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FIG. 7: The nonlinearity parameter /nl versus fio-dcc and pin: numerical results and sudden decay approximation at first and 
second order. 



B. Numerical solutions 

The transfer parameter at first order relating the initial curvature perturbation on uniform curvaton density hy- 
persurfaces to the final value of the total curvature perturbation is defined as [3, S 



ri 



Cl,fin 



(5.10) 



We define the transfer parameter at second order 



C2,fin 



(5.11) 



relating the final value of the total curvature perturbation to the initial curvature perturbation on uniform curvaton 
slices. The values for ri and r2 coincide for our choice of initial condition, Cio-jin = 1, with the final values of the 
curvature perturbations at first and second order, Ci.fln and C2,fin, and are given in Fig. versus rio-doc and pi^. 

The non-linearitjAparameter using first order perturbation theory is given in terms of the transfer parameter defined 
inEq. (|5lT)|l as (iQ 



/nl - -j — . 



(5.12) 



Using second order perturbation theory we find the non-linearity parameter /nl from Eq. H5.1|l . expanding C, to 
second order, and get in terms of the transfer parameters at first and second order 



5 5 r2 5 
4ri Gri^ 3 



/nl = -jz~ + ^ TT ~ o 



(5.13) 



In the above calculations we identified <^g with the part of Ci linear in the curvaton field fluctuation, i.e. the first term 
in Eq. (jO|) . 

We can finally relate the transfer parameters ri and r2 to the total curvature perturbation, (^ = Ci + 2^2: evaluated 
after the curvaton has decayed, 



Cfin = riClaM + -jr2Cla.} 



(5.14) 



Using the definition of the curvature perturbation on uniform curvaton density hypersurfaces, Eq. H2.19|l . and the 
expression for the curvaton energy density in terms of the curvaton amplitude, Eq. (|5.4(l . we get 



^ 2 5a. I ( 2 \ fSai 

3 o-Q 3 V 3 / V (To 



(5.15) 
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FIG. 8: The nonlinearity parameter /nl versus ficrdec and pin: numerical results and sudden decay approximation at first and 
second order, detail of Fig. 



We can compare these results to the expression found using the AiV-formalism, expressed in terms of the curvaton 
perturbation (instead of the normalised energy perturbation used in Section [lV|l |l2l.l88l|. 



CsB = N^^Scr + -N^^^Sa^ . 
The curvature perturbation employed in the A A^- formalism, C2SB , is related to ( by Eq. 12.32|l , and we get 



2 5ai 1 
CsB = 0^1 ^ o I ''I + q''2 - 0*^1 

d (To O V -J O 






and therefore 
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(5.16) 



(5.17) 



(5.18) 



C. Curvaton amplitude evolution 



So far we assumed that the curvaton field doesn't evolve between the end of inflation and curvaton decay. In order 
to allow for the evolution of curvaton field amplitude a we assume that it depends on the initial value set during 
inflation CTc by ct = cr((To) which gives for the curvaton field fluctuation |Mll0Llia| 



1 



5a = a 5(7 o H — cr 5a, 
2 



"iC_2 
e t 



where a' = da /dac- For the first order energy density we then find (including again the quadratic term) 

2 / // \ 



POcr ctq 

We therefore get for the non-linearity parameter 



5aci 



1 + c^o— 2 o-Q 



/2 
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(5.19) 



(5.20) 
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NL 
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l + CTo- 
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(5.21) 



instead of Eq. (|5.13() above. However, in order to calculate a numerical value for /nl we now have to calculate the 
evolution of a in detail and specify a curvaton model. We shall therefore not pursue this issue further and refer to 
Ref. llOl where this issue was studied in detail. 
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FIG. 9: The difference of tlie numerical and the sudden decay approximation value of the non-linearity parameter, A/nl, 
versus ficdcc and pin. 



D. Results and discussion 

Our results are summed up in Figs. l7li?l in Figs.E|and|Slwe plot the non-linearity parameter calculated numerically 
using Eq. H5.12|l at first order and Eq. (|5.13l) at second order for two different parametrisations, namely fio-dccand 
Pin- In the same figures we also plot the non-linearity parameter /nl in the sudden decay approximation at first and 
second order from Eqs. (|5.8|l and H5.9|l . respectively. We truncated the graphs at /nl = 120 in accordance with the 
current observational bounds (see below). We first note how well the sudden decay approximation and the numerical 
solution at first and second order agree. However, at first order we find that /nl = 5/4 for ^adcc = 1 or large values 
of Pin, whereas including the second order effects we get /nl = —5/4 for Qa-dcc ~ 1 or large values of pin. 

In Fig. El we plot the difference between the non-linearity parameter obtained numerically and using the sudden 
decay approximation, both at second order. 



A/nl = /] 



NL 



numerical 



/] 



NL 



suddcndccay 



(5.22) 



We now see more clearly the excellent agreement of the sudden decay approximation for large parameters pin and 
fio-dec- However, for small pin and ilo-doc the sudden decay approximation works less well, deviating from the numerical 
solution by up to 10%. 

We finally give the observational constraint on the non-linearity parameter /nl from the recently published WMAP 
three-year data. Spergel et al. found |43: —54 < /nl < 114 (at 95% confidence level). The curvaton model is 
therefore well within the current observational bounds. However, if future observations give a large negative value for 
the non-linearity parameter, the curvaton model would be ruled out, at least without strong evolution of the curvaton 
amplitude from the beginning of the oscillations to curvaton decay, as pointed in in Section IV CI above. 

Acknowledgments 

The authors are grateful to Jussi Valiviita and David Wands for useful comments. KAM is supported by PPARC 
grant PPA/G/S/2002/00098, DHL is supported by PPARC grants PPA/V/S/2003/00104, PPA/G/O/2002/00098 
and PPA/S/2002/00272 and EU grant MRTN-CT-2004-503369. Algebraic computations of tensor components were 
performed using the GRTensorII package for Maple. 

APPENDIX A: GOVERNING EQUATIONS 

Here we first give the governing equations on large scales in the general case without gauge restrictions and then 
the equations given in Section IlIII in terms of non-normalised energy densities. 



1. Governing equations wthout gauge restriction 



In this subsection we give the governing equations on large scales in the general case without any gauge restrictions, 
i.e. without choosing a particular hypersurface. 



a. First order 
Energy conservation of the a-fluid is given from Eq. H2.4I) at first order as 

Spia + iH {5pia + SPia) - ilpl {pOa + Poa) " Qoq^I " SQia = . (Al) 

Total energy conservation follows from Eq. (|A1|I above, and using Eqs. H2.5() and H2.16|l . is given by 

6pi + 3H {5pi +SPi)~ SV-i (po + Po)^0. (A2) 

The — Einstein equation is given from Eq. (|2.1|) as 

3H {H(J)i + VJi) = -^TiGSpi . (A3) 

b. Second order 
Energy conservation of the a-fluid is given from Eq. p. 4(1 at second order as 

dp2a + 3if {Sp2a + SP2a) - 3?/'2 {pOa + Poa) - 6 {Spia + 6Pia) i'l - 12 (poa + Poa) V'lV'l 

- Qoc (02 - 0?) - 20i'5Qic - (JQaa = , (A4) 

and, following a similar route as at first order, the conservation of the total energy density is given at second order by 
Sp2 + 3H {5p2 + 5P2) - 3^2 {po + Po)-Q {Spi + SPi) ^1 - 12 (po + Po) Mi = , (A5) 

and the — Einstein equation is given by 

3H^ (02 - m) + 3 (iJ^2 - i^l) + l2H^i (V'l - 0i) = -^ttG5p2 ■ (A6) 



2. Governing equations in terms of non- normalised energy densities 

In this subsection we give the governing equations presented in Sections IIII Bl and IIII CI above in terms of the 
normalised quantities in terms of the non- normalised energy densities and decay rate. We use the number of e-foldings 
N as time coordinate and work throughout in the flat gauge (omitting the "tilde"). 



We get at first order 



and at second order 



Vi.+ (3 + §)<5pi.-i§^^Pi=0, (A7) 

Vi^ + 4<5pi^ - i-<5pi. + \^^5p^ = , (A8) 

n Z H Pq 



Sp'2.+ (^ + ^)sp2.-^'-^Sp,. + ^^(l'-^-bp2)^0, (A9) 

\ HJ H Po H pq \A Po 2 J 

5p'2,^^Sp2, - ^5p2. + ^'-^Sp,^ - ^^ (I'J^ - lsp2) =0. (AlO) 

' H H Po H Po \4: Po 2 J 
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